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Abstract The ability to perform quantum error correction is a significant hurdle 
for scalable quantum information processing. A key requirement for multiple-round 
quantum error correction is the ability to dynamically extract entropy from ancilla 
qubits. Heat-bath algorithmic cooling is a method that uses quantum logic operations 
to move entropy from one subsystem to another, and permits cooling of a spin qubit 
below the closed system (Shannon) bound. Gamma-irradiated, ^^C-labeled malonic 
acid provides up to 5 spin qubits: 1 spin-half electron and 4 spin-half nuclei. The nu¬ 
clei are strongly hyperfine coupled to the electron and can be controlled either by ex¬ 
ploiting the anisotropic part of the hyperfine interaction or by using pulsed electron- 
nuclear double resonance (ENDOR) techniques. The electron connects the nuclei to 
a heat-bath with a much colder effective temperature determined by the electron’s 
thermal spin polarization. By accurately determining the full spin Hamiltonian and 
performing realistic algorithmic simulations, we show that an experimental demon¬ 
stration of heat-bath algorithmic cooling beyond the Shannon bound is feasible in 
both 3-qubit and 5-qubit variants of this spin system. Similar techniques could be 
useful for polarizing nuclei in molecular or crystalline systems that allow for non¬ 
equilibrium optical polarization of the electron spin. 
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1 Introduction 

Quantum Error Correction (QEC) is a critical tool for protecting quantum informa¬ 
tion against the imperfections of realistic devices and to scale quantum processors up 
to many qubits. Although a well developed theory exists mu and experimental real¬ 
izations at the several qubit level are recently emerging Bia , there remain challenges 
for many potential implementations. Nuclear Magnetic Resonance (NMR) quantum 
information processing has demonstrated a high degree of quantum control mm 
and the ability to efficiently characterize the noise, both intrinsic and extrinsic, that 
affects the fidelity of a quantum process ifT^ . One piece has however been missing: 
the ability to efficiently polarize nuclear spin qubits on demand. The threshold theo¬ 
rem for quantum computation tells us that a quantum circuit can be simulated with 
arbitrarily high precision using a polynomial amount of resources as long as the error 
per gate p is below a certain threshold value pth is . The theorem relies on assuming 
that the ancilla qubits are in pure states at the beginning of each cycle of fault-tolerant 
QEC. Eor example, the first layer of concatenation of QEC typically reduces the ef¬ 
fective error rate from p to cp^, where c is a system-dependent constant. However, 
this theoretical gain is not generally achieved for the impure ancilla qubits charac¬ 
teristic of real implementations im. Thus, an efficient and experimentally feasible 
method for cooling qubits to high purity prior to each QEC cycle is desirable for all 
circuit implementations, including those based on nuclear spins. 

Heat bath algorithmic cooling (HBAC) |[T3llT5l[T7lfT^ is an efficient method for 
extracting entropy from a set of system qubits, allowing qubits to be cooled below 
the bath temperature (i.e. beyond the closed-system, or Shannon, bound). Solid state 
NMR experiments have demonstrated a sufficient level of coherent control to execute 
multiple rounds of algorithmic cooling, leading to spin polarizations exceeding the 
thermal polarization |[T^l20ll . However, in a typical NMR setup, very low spin polar¬ 
ization at thermal equilibrium will require highly precise control over tens of nuclear 
spin qubits in order to polarize one ancilla qubit to order unity. Here, we explore the 
use of an electron spin to assist in the HBAC protocol. Due to its much larger gyro- 
magnetic ratio compared to nuclei, the electronic thermal spin polarization is about 
10^ times larger, and spin-lattice relaxation rates scale by a similar factor. Exploiting 
the electronic spin-lattice relaxation as a reset operation, the electron can connect a 
set of nuclear spins to heat bath with an effective temperature much lower than the 
equilibrium nuclear spin temperature. 

The stable malonyl radical CH(COOH )2 has been extensively studied by elec¬ 
tron spin resonance (ESR) studies of gamma-irradiated single crystals. The hyperfine 
tensors of the a-proton and of the ^^C-labeled methylene carbon were previously 
published ll2TH2^ . However, the hyperfine tensors of the ^^C-labeled carboxyl car¬ 
bons have not yet been reported. Including these two carboxyl carbons and assuming 
they are spectroscopically distinct, the molecule can in principle realize a 5-qubit en¬ 
semble quantum information processor, with 1 electron and 4 nuclear spins. In this 
paper we study the per-^^C-labeled radical and determine the hyperfine tensors of 
the carboxyl carbons via ESR and electron-nuclear double resonance (ENDOR) ex¬ 
periments on single crystal samples at room temperature. The carboxyl tensors are 
in fact different owing to different dihedral angles of COOH relative to the C-C- 
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C plane to accommodate the hydrogen bonding network; the carboxyl group with 
a larger dihedral angle has a slightly weaker carbon hyperfine coupling. The previ¬ 
ously published tensors describing the electronic g-factor and proton and methylene 
carbon hyperfine couplings are confirmed by our measurements. Given full knowl¬ 
edge of the spin Hamiltonian, we determine the optimal magnetic field orientations 
for carrying out quantum algorithms, focusing on HBAC as the example of inter¬ 
est. We consider two distinct methods for obtaining high fidelity coherent control, (i) 
use of anisotropic hyperfine interactions ll24H26ll facilitated by GRAPE Ell pulses, 
and (ii) pulsed ENDOR control sequences. Eor brevity, we refer to approach (i) as 
Anisotropic Hyperfine Control (AHC). The two control schemes have very different 
criteria for choosing a suitable magnetic field orientation, given the desire to reduce 
the durations and optimize the fidelities of the quantum operations. Simulations of 
HBAC with realistic values of electronic spin-lattice relaxation time (Ti), intrinsic 
dephasing time (T 2 ) and ensemble dephasing time (T^) are carried out for both the 
3-qubit (methylene and 5-qubit (per-^^C labeled) versions. The results indicate 
that experimental cooling of nuclear spins beyond the thermal electron spin temper¬ 
ature, while challenging, is feasible in this quantum processor. These techniques can 
be applied to malonic acid or similar radical molecules at very low temperatures, or to 
optically polarizable systems such as NV centers in diamond ll28H30ll or molecules 
with photoexcited triplet states ED at room temperature, to reach the nuclear spin 
polarizations of order unity useful for QEC. 

The reminder of the paper is organized as follows. Section provides relevant 
information about the properties and preparation of the single crystal samples. In 
section we review the extraction of g-factor and hyperfine tensors from orienta¬ 
tion studies of ESR and ENDOR transitions, and report the tensors we determined 
including the carboxyl tensors. Section [^presents criteria for choosing optimal 
orientations of the magnetic field with respect to the crystalline axes, to facilitate 
control using either anisotropic hyperfine or pulsed ENDOR techniques. Realistic 
simulations of HBAC protocols for both 3 and 5 qubits are described in section 
Conclusions are drawn in section | 6 l 


2 Sample preparation 

Malonic acid, CH 2 (COOH) 2 , crystalizes with a triclinic unit cell and belongs to the 
PI space group |[32j[33l at temperatures above 47 K. At lower temperatures, a struc¬ 
tural phase transition occurs that has been discussed previously 13414361 . Above 47 
K, there are two molecules per unit cell related by inversion symmetry, making them 
magnetically equivalent. We denote the methylene and carboxylic carbons as Cm 
and Cl^ 2 , respectively. A schematic of the radical, obtained by removing one of the 
methylene protons, is shown in Eig. 1. The unpaired electron is of p-orbital character 
|22l [37l[38l . Malonic acid powder with all possible isotopic labelling configura¬ 
tions was purchased either from Sigma Aldrich or Cambridge Isotopes. Single crys¬ 
tals were grown by slow evaporation from aqueous solutions at room temperature. To 
form radicals, crystals were irradiated to a dose of about 2 kCy at room temperature 
with 7 -rays from a cobalt-60 pencil. Annealing at 60° for 12-15 hours following the 
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irradiation suppresses ESR signals from all other radical species except for the most 
stable radical, CH(COOH )2 EtI. Depending on the labelling configuration, 2, 3, 
4 or 5 qubit samples are obtained as (e-H), (e-H-C^), (e-H-Ci^2) and (e-H-C^-Ci^2) 
respectively, where e denotes the electron spin. The j3-protons of the carboxyl groups 
contribute to ESR line broadening and in some orientations can give rise to observ¬ 
able splittings 13^ . However, due to weak hyperfine coupling (< 6 MHz isotropic 
coupling |[39l ) they are not useful as qubits. 




Fig. 1 (a) Molecular structure of CH(COOH )2 with unpaired electron density distribution schematically 
represented by the blue shaded region, x and z are two principal axes of the a-proton hyperfine tensor, with 
z along the Cm-Of-proton interatomic vector and x along the cylindrical symmetry axis of the electronic 
p-orbital. The y axis, not shown, is nearly parallel to the C 1 -C 2 interatomic vector. The direction of the 
static field Bq can be described using polar angle 6 and azimuthal angle (j) in the principal axis system of 
the a-proton hyperfine tensor. The directions of the crystallographic axes a, b and c with respect to the 
molecular structure are illustrated by the green axes (^IdOj. c lies very close to the y axis of the a-proton 
principal axis system. The direction cosines of a, b and c in the a-proton principal axis system are (0.1426, 
-0.6588, 0.7387), (-0.9729, -0.1888, 0.1335) and (0.0222, 0.9969, -0.0752), respectively, (b) Schematic of 
the crystal mounted on a sample holder for the orientation study (see sectionPrimed axes (x', y' and 
z') indicate the lab frame in which measurements are taken. The curved green arrow indicates rotation of 
the crystal for orientation studies. 



3 Spin Hamiltonian 

The spin Hamiltonian of malonic acid with K nuclei contains electron Zeeman, electron- 
nuclear hyperfine, and nuclear Zeeman terms as following: 

= Jfe + ^hf + ^n= I^BgapBoaSp + ^ [iTtAlpSalp - Ynl^Boaj . (1) 

n=l 

Here a, j3 G {v,y,z}, and the repeated index implies summation over all the values of 
the index (similar to the Einstein summation convention, except all indices appear as 
lower indices). We scth = l so that all Hamiltonians will appear in angular frequency 
units. S(^ p and ^ represent electron and nuclear spin operators, is Bohr mag¬ 
neton, Bq = {Box^Boy^Boz) is external magnetic field and is the gyromagnetic ratio 
for nuclear spin n. The second rank tensors g whose elements are and whose 
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elements are , are the electron g-tensor and the hyperfine tensor describing cou¬ 
pling to nuclear spin n. In X-band (^10 GHz) ESR, the electron Zeeman interaction 
is the dominant Hamiltonian term. The nuclear dipole-dipole interaction is neglected 
since it is typically at least two orders of magnitude smaller than the hyperfine inter¬ 
action. When the nuclear Zeeman and hyperfine interaction energies are comparable 
and much smaller than the electron Zeeman energy, taking Bq = the secular spin 
Hamiltonian can be written: 

^ = (OsS, - f + 2kS, f + 5"/^). (2) 

n=\ n=\ 

Here, CDs = SzzI^bBq is the electron Larmor frequency, cOj = Yn^o is the nuclear Lar- 
mor frequency, and = yj{A^Y + 


3.1 Tensor extraction from ESR and ENDOR data 


In this section, we briefly review the Hamiltonian determination method explained 
in ED. Eirst, we will consider the procedure to determine the g-tensor. Eor a static 
magnetic field Bq = lyih) QxJy and 4 are direction cosines of the magnetic field 
in the axes of the crystal, and -h /^ = 1), the electron Zeeman energy levels can 
be rewritten as 

E± = ±^-HBBo{laraplpy^\ (3) 

where we introduce f = g • g. In the absence of nuclear spins, the magnetic resonance 
condition AE = hv = leads to 


hv 

I^bBq 


AE _E^-E- 
I^bBq JD^bBq 




(4) 


where g is the experimentally observed value, AE" is the difference between and 
E-, and V is the frequency of the microwave (mw) field. Eor Bq = (cos 0, sin 0,0), 
using the symmetry the dependence of g on 0 can be written as 

= la^af^h = 0 ^ Eyy sin^ 0 COS 0 sinO . (5) 

By measuring g(0) at various 0 and solving Eq. ^,rjcx,ryy and E^y can be deter¬ 
mined. Similarly, by orienting Bq within the v-z and y-z planes and measuring g{6) 
for various 6 (where 6 is the angle between B and z, or between B and y), and 

Fyz can be determined. Therefore full knowledge of V can be obtained. 

Now we consider the electron-nuclear coupled spin Hamiltonian which contains 
and In this case, the 1st order approximation of the energy levels of the 
system are 


K A ^2 y2 d2 

EmsM, =8BbBoMs + h ) + 

n=i S lyis 

- ^r«5o(/a(g- 


( 6 ) 
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where = g • A • A • g, and Mg and are quantum numbers for the electron spin and 
the nuclear spin, respectively. All ESR transitions can be classified into two types: 
allowed and forbidden transitions. Allowed transitions correspond to electron-only 
spin flips, while forbidden transitions correspond to the flips of electron spin along 
with one or more nuclear spins. In the absence of anisotropic hyperfine coupling (the 
terms with coefficient in Eq. Q), forbidden transitions are completely suppressed 
and produce no observable ESR signal. When ^ 0, these transitions can lead to 
signals comparable to, but typically smaller than, the allowed transitions. Eor the pur¬ 
pose of hyperfine tensor determination, it suffices to track only allowed ESR transi¬ 
tions. The nuclear Zeeman energy is generally much smaller than hyperfine coupling 
strength and is usually omitted in the traditional process of determining hyperfine 
tensors using ESR spectra. However when we determined hyperfine tensors in this 
work, after the traditional process we implemented an optimization process which 
takes the nuclear Zeeman energy into consideration. We will discuss this further in 
Section 3.2 In the present section we will stick to the traditional process. Therefore, 
after omitting the nuclear Zeeman energy, the ESR allowed transition energy gap in 
the 1 electron-1 nucleus system is 


AE = s^bBo + 


(7) 


where the nuclear spin quantum numbers M/ = ±1/2 represent spin up and spin down 
states, respectively. The ESR spectrum then splits into two distinct peaks separated 
by a frequency A = Similar to Eq. (j^, when the static magnetic 

field is Bo = (cos 0, sin 0,0) we have 

(g(0)A(0))^ = cos^ O-hTy^ sin^ 0 ± cos 0 sin 0, (8) 


and, as before, expressions for Bq in the v-z and y-z planes relate the other tensor 
components of to (gA)^. 

Another method to determine hyperfine tensors is to make use of ENDOR exper¬ 
iments. Eor ENDOR transitions, which correspond to AMg = 0 and AM/ = ±1, the 
energy gap in the 1 electron-1 nucleus system can also be readily obtained from Eq. 
•HI* and its square is given by 


AE^ 


= fBl 




(lar^ph)- 


AnMgjBQ 

g 




(9) 


Therefore, the ENDOR spectral peaks also split into two for Mg = ±.\/2. Supposing 
that the observed ENDOR frequencies are v± Fox Mg = ±1/2, Eq. (|^ can be rewritten 
as 

= ^/ + ^ (^«±3 ) =F y (^« (g • A) ajS ^(3), (10) 

where V/ = yB^l2%. An equation similar to Eqs. and ^ can be obtained when the 
static magnetic field is Bq = (cos 0, sin 0,0), 

±(v^ -v+) = (g-A)^;(Cos^e + (g-A) 3 ,j,sin^e + 2(g-A);cj,cosesine. (11) 

ZV[ 

Therefore, the magnitude of A can be extracted by measuring the orientation depen¬ 
dence of g(v?: — v^)/2v/ in three distinct rotation planes. 
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3.2 Continuous-wave ESR results 

The tensor extraction model relies on two assumptions: the three planes of measure¬ 
ment are mutually non-parallel and the axis of the rotation to go from one to another 
belongs to both planes. The orientation experiments were designed to ensure that 
both assumptions are satisfied. In each of the three non-parallel planes, we took 24, 
9 and 46 measurements with 8°, 20° and 4° angle steps for the methylene-labeled, 
carboxyl-labeled and fully-labeled MA samples, respectively. The methylene-labeled 
data was used to extract the a-^H and tensors; the carboxyl-labeled data was 
used to extract the (average) tensor, and the fully-labeled data was used to 

confirm all the tensors by fitting with simulated spectra. As discussed previously in 
Section |3.1| we need only to detect the spectral positions of the allowed transitions, 
which can be distinguished from forbidden transition peaks because of their larger 
intensities. Each measured spectrum is fit to a set of allowed transition peaks defined 
by their amplitude, frequency and a common line width. 

The peak positions for all measurements in one plane give a set of trajectories for 
that plane, from which we obtained the dependence of g and on 0. Doing this in 
all three planes allows extraction of the g-factor tensor g and hyperfine tensors 
by solving Eqs. ^ and After this, we used the method of least square fitting 
to optimize g and by minimizing the difference between the experimental peak 
trajectories and the simulated trajectories. It should be noted that in X-band ESR the 
nuclear Zeeman energy for is about 15 MHz, which is comparable to the hyperfine 
coupling of a-^H in some orientations. Therefore, in order to obtain more accurate 
A'^ tensors, we took into account the nuclear Zeeman energy when generating the 
simulated peak trajectories in the optimization procedure. The only noticeable differ¬ 
ence from including this energy was for the a-^H, as expected. 

The g, A*^, A^'^ and A^i ’2 tensors were determined from continuous (CW) ESR 
spectra as described above, and are listed in TableThe hyperfine tensors of Ci and 
C 2 are similar and require ENDOR measurements to be distinguished (see next sec¬ 
tion). Erom the ESR data we obtained an average hyperfine tensor describing Ci and 
C 2 , denoted as A^i ’2 Negative signs of the principal values for A^i ’2 determined 
by a quantum chemical calculation using the Gaussian program 1461 . Similar to the 
case of a-^H ||22l [371 |38l, the negative principal values are due to a negative spin 
density on Ci and C 2 . 

The results for g, A*^ and A^'^ are consistent with the published results in 1271 - 
[231 EH [3711401. The simulated ESR peak trajectories generated from the tensors in 
Tabled are in excellent agreement with the experimental data, as the examples in 
Eigs. [^nd[^ demonstrate. 


3.3 Continuous-wave ENDOR results 


In this section we describe the use of ENDOR, which has higher spectral resolution 
than ESR, to extract the distinct tensors describing A^i and A ^2 From Eq. (10) we 
see that for each nucleus there are two resonant frequencies v_ and v+. However, 
in experiments we found that one or sometimes two of the four peaks of Ci and C 2 
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Direction cosines to A*^ principal axis 



X 

Y 

Z 

gxx 

2.00250 ± 0.00038 

-0.1657 

0.9779 

0.1272 

gyy 

2.00373 ± 0.00037 

-0.9811 

-0.1766 

0.0797 

gzz 

2.00417 ± 0.00036 

0.1004 

-0.1115 

0.9887 

A» 

-26.6 it 2.8 

1 

0 

0 

AH 

-56.0 it 0.7 

0 

1 

0 

^zz 

-91.5 ±0.6 

0 

0 

1 

A Cm 
^xx 

24.5 ± 1.0 

0.0696 

-0.0019 

0.9976 

A Cm 

^yy 

43.0 ± 1.3 

0.9962 

0.0530 

-0.0694 

A Cm 
^zz 

212.3 ±0.6 

-0.0528 

0.9986 

0.0056 

aCi,2 

-^XX 

-36.1 ±0.3 

0.0627 

0.0082 

0.9980 

A^1^2 

n.yy 

-39.3 ± 0.3 

0.9805 

-0.1870 

-0.0601 

aCi,2 

^ZZ 

-40.6 ± 0.3 

0.1862 

0.9823 

-0.0198 


Table 1 Electronic g-factor and hyperfine coupling tensors determined from CW ESR measurements 
(principal hyperfine values are given in MHz). Principal values are given in the left column, while direction 
cosines relative to the principle axis system of the a-proton are given in the three right columns. The 
uncertainties reflect a 90% confidence interval. The tensor gives an estimate of the average of 
and A ^2 



Fig. 2 Comparison between the simulated ESR peak trajectories (red lines) and the experimental data 
(blue dots) of the methylene ^^C-labeled sample in one of the planes measured. The direction cosines of 
the normal of this plane are (-0.3376 0.9405 0.0391) in the principal axis system of A*^. 


were obscured by the large peaks around 14.5 MHz that come from spins which 
have small hyperfine couplings with the electron radical (Fig. [^. The two higher 
frequency peaks (in the range of 20 — 25 MHz) were clearly resolved in all measured 
orientations. Judging from the sign of A^T 2 ^ we infer that these two peaks correspond 
to the and frequencies, with 
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Fig. 3 Comparison between the simulated ESR peak trajectories (red lines) and the experimental data 
(blue dots) of the per-^^C-labeled sample in one of the planes measured. The direction cosines of the 
normal of this plane are (-0.3678 0.9287 -0.0483) in the principal axis system of 


where n= 1,2 for Ci and C 2 and — Vjgg • Using only the v+ fre¬ 

quencies requires a different approach than directly using Eq. O- Instead of mea¬ 
suring the orientation dependence of we measured that of 

^2((y^)2 _ (yn)2) in thi-ee non-parallel planes to extract using Eq. (12), and then 
used to calculate and A ^2 each plane, we took 9 measurements with 20° 
angle steps for the carboxyl-labeled MA sample. It should be mentioned that by using 
Eq. (12) to obtain A^i and A^ 2 ^ neglected the anisotropy of g. Eor the hyperfine 
tensors, this approximation causes a difference in the principal values less than the 
experimental uncertainty of ^ 0.1 MHz. The experimentally determined A^i and A ^2 
tensors are listed in Table O 




Direction cosines to A*^ principal axis 



X 

Y 

Z 


-37.0 ±0.1 

-0.5310 

-0.0010 

0.8474 

/\yy 

-40.5 ±0.1 

0.8452 

-0.0724 

0.5295 

A^i 

^zz 

-43.6 ±0.1 

0.0608 

0.9974 

0.0392 

-^xx 

-34.0 ±0.1 

-0.3477 

-0.0449 

-0.9365 

A^2 

/\yy 

-37.4 ±0.1 

-0.2478 

0.9677 

0.0456 

A^2 

^zz 

-39.6 ±0.1 

-0.9043 

-0.2479 

0.3476 


Table 2 Hyperfine coupling tensors for and A ^2 determined from CW ENDOR measurements (prin¬ 
cipal hyperfine values are given in MHz). Principal values are given in the left column, while direction 
cosines relative to the principle axis system of the alpha-proton are given in the three right columns. The 
uncertainties reflect a 90% confidence interval. 


Figs. 0 and show the experimental ENDOR and ESR spectra, together with 
the simulated spectra generated using the tensors given in Table for the carboxyl- 
labeled sample when the orientation of the magnetic field is (0.6156 -0.7179 -0.3249) 
in the principal axis system of AIn Eig.|^ we can see that in this orientation there 
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are three observable peaks of Ci and C 2 , while the peak with the lowest frequency 
is obscured by the peaks around 14.5 MHz. Excellent agreement between experi¬ 
ment and simulation can be seen in these examples, and a similar level of agreement 
was found for spectra in all measured orientations, indicating that the tensors in Table 
[ 2 ] are accurate. 



Fig. 4 Comparison between the simulated ENDOR spectrum (red line) and the experimental spectrum 
(black line) of the carboxyl ^^C-labeled sample. The orientation of magnetic field is (0.6156, -0.7179, - 
0.3249) in the principal axis system of . In the area around 14.5 MHz, there are large peaks which come 
from the distant spins that have small hyperfine couplings with the electron radical. The ENDOR peak 
for C 2 with the frequency = 14.9 MHz is obscured by these ^H peaks. The other three ENDOR peaks 
for Cl and C 2 with frequencies v]_ = 24.1 MHz, = 22.1 MHz and vi = 16.8 MHz are clear. 



Fig. 5 Comparison between the simulated ESR spectrum (red line) and the experimental spectrum (black 
line) of the carboxyl ^^C-labeled sample. The orientation of magnetic field is (0.6156, -0.7179, -0.3249) in 
the principal axis system of . Using the tensors given in Table|^ the ESR peak positions are reproduced 
well and the line broadening due to the hyperfine couplings of the two carboxyl carbons fits well in all 
measured orientations. 
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4 Optimizing magnetic field orientation for quantum control 

4.1 Orientation criteria for AHC 

In the presence of Bq = Bqz, a hyperfine-coupled nuclear spin is quantized along the 
direction of an effective field 

Beff = (5o±:7rA,,/7)z±(;r5/7)x, (13) 

where B = and the ± sign depends on whether the electron spin is 

parallel (|t)) or anti-parallel (||)) to Bq. The angle of nuclear spin quantization 
axis from z is denoted as = arctan {—B^/-\-C0j/7r)) for electron spin |t) and 
7]^ = arctan {—B^/— a>j jn)) for electron spin ||). 

A requirement for achieving high control fidelity is to implement gates on a 
time scale fast compared to the electronic dephasing time T^. The electron spin con¬ 
trol Hamiltonian in the rotating frame is = mi(r)(5xcos(0(r)) Sy sin 

where coi and (j) represent the microwave amplitude and phase, respectively. Consider 
(Di = constant, and 0 = 0. Nuclear spin flips can occur via electron-nuclear flip-flip 
forbidden transitions driven at a rate C0ism(ri^)/2 where rj^ = (r]^ — T]p/2. The 
on-resonance allowed ESR transition is driven at tOi cos(r]'^)/2. Depending on the 
nuclear isotope and the dc field orientation, sin(r ]'^)/2 can be very small, preventing 
efficient control of the nuclear spins via microwave excitation. Thus, we aim to find 
an orientation that yields the set {A^^B^} that maximizes the ratio of forbidden to 
allowed transition rates, tan(r 7 '^). Secondly, all allowed and forbidden transition fre¬ 
quencies should be separated from each other by at least the ESR line width (^12 
MHz in the malonyl radical) in order to achieve universal control. 

Eor the optimal orientation search, is fixed at 2;r • 10 GHz/y^ « 3568 G (note 
7e = gjiB) and 6 and 0 (see Eig. [^a)) are varied at 1° angle steps. The maximum 
values of tan(r 7 '^) are 1.3, 0.14, 0.016, and 0.017 for H, C^, Ci, and C 2 , respec¬ 
tively. The duration of a GRAPE optimal control pulse to accomplish an arbitrary 
unitary operation will scale as l/tan(r]'^). We have measured a nearly temperature- 
independent electron spin dephasing time ^ 5 /is in this system, due to hyperfine 
coupling with distant protons, and their mutual dipolar interaction induced dynamics. 
Given that the shortest GRAPE pulses for arbitrary gates involving H require ^500 
ns, it is clear that high fidelity control of Ci, and C 2 via AHC is simply not possi¬ 
ble in this system. Moreover, the 5-qubit system contains 80 ESR resonances over 
a maximum spectral range of 250 MHz, consisting of 16 allowed and 64 forbidden 
transitions. Consequently, we find max^ 0 {min( 5 )} = 1.4 MHz, where the elements 
of 5 are the distances between any two ESR transitions and max^ indicates the 
maximization over all sets of {0,0} with 1° angle step. In practice, it only makes 
sense therefore to focus on the 3-qubit sample (e-H-C^) for optimizing AHC control. 
Nonetheless, for completeness we show in Table the best possible orientation for 
the 5-qubit system given these criteria. 

Eor the 3-qubit sample, there are 4 allowed transitions and 8 forbidden transitions. 
It is still not possible to find orientations where min(5) >12 MHz, rather we find that 
maxe 0 {min( 5 )} = 7.6 MHz. We focus on orientations that satisfy min(5) > 6 MHz, 
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and then choose an orientation in which tan(r7^) and tan(r7^"') are both as large as 
possible. The forbidden transition rates at this orientation are 0.1418 and 0.1203 for 
H and C^, respectively, and min(5) =6.1 MHz. These results are summarized in the 
table below: 




A’' 

?^(MHz), 5^(MHz), tan(77^) 



0,0 

H 

C™ 

Cl 

C 2 

min(5) 

(MHz) 



-62.56, 

121.23, 

-40.03, 

-37.29, 


5 qubit 

54°,123° 

23.67, 

90.54, 

3.27, 

1.63, 

1.1 



0.1948 

0.0303 

0.0161 

0.0093 




-76.60, 

29.21, 




3 qubit 

29°, 13° 

27.07, 

17.78, 



6.1 



0.1418 

0.1203 





Table 3 Optimal orientations for AHC experiment with the fully-labeled (5 qubit) and -labeled (3 
qubit) malonyl radical, and corresponding values of hyperfine coupling constants and 5”, and the quan¬ 
tity tan(r7”) characterizing the controllability of nuclear spin n. min(5) represents the minimum separation 
between ESR transitions. 


4.2 Orientation criteria for ENDOR 


In the ENDOR control scheme, nuclear spin rotations are implemented directly by ap¬ 
plying radio-frequency (RE) pulses on resonance with the nuclear transitions. Hence, 
tan(T7) does not need to be large. In fact, it is desirable to minimize tan(r]), because 
the electron spin-lattice relaxation process induces nuclear spin relaxation through 
the two-spin flip forbidden transitions. As before, the allowed ESR transitions must 
be separated by at least the ESR line width, and the NMR transitions must be sep¬ 
arated by at least the NMR (ENDOR) linewidth. The latter is determined by the 
nuclear dephasing time and is therefore much narrower than the ESR 

line width. The disadvantage of ENDOR control is that the RE pulses are typically 
of order yehn ^ 10^ times slower than the microwave control of electron spin. This 
severely limits the ability to perform arbitrary quantum algorithms within the elec¬ 
tron however, in algorithmic cooling the electron spin is always in an eigenstate 
during the nuclear rotations so that electronic dephasing is not an issue. Since the 
hyperfine tensors of Ci and C 2 are very similar, they cannot be separately addressed 
by microwave pulses, but only with RE pulses. At first glance, this seems to forbid 
selective swap gates between the electron and Ci or C 2 , which is an essential step 


for HBAC. Eortunately, we show in Section [5^ that it is possible to realize a pulse 
sequence combined with the electron refresh step in order to polarize both Ci and C 2 
to the bath polarization, as long as all NMR transitions can be addressed selectively, 
and the nuclei H and Cm are separately addressable by microwave pulses. 

The RE pulse durations are typically on the order of tens of /is, corresponding 
to excitation bandwidths less than 100 kHz. Thus, for the 5 qubit sample, we search 
for an orientation in which all NMR transition frequencies are at least 500 kHz apart. 
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and all allowed ESR transition frequencies with the exception of the Ci and C 2 res¬ 
onances are at least 12 MHz apart. Then among this set, we select one orientation 
in which tan(r]) is as small as possible for all nuclei. The optimal orientation for the 
fully-labeled molecule and relevant parameters in that orientation is given in Table 
m The minimum distance between ESR transitions in this orientation is 6 MHz due 
to Cl and C 2 . All other ESR transitions are at least 13 MHz apart from each other. 
Some forbidden transitions may be very close to the ESR allowed frequencies, but 
this is not a serious issue since the forbidden transitions can be ignored when tan(r]) 
is small. 

Eor the 3 qubit molecule, it is relatively easy to find orientations in which the al¬ 
lowed ESR transitions are separated by at least 12 MHz. Among that set, an optimal 
orientation is chosen for which tan(r 7 ) is minimized for both H and C^. The results 
are summarized in Table [U 





f^(MHz), 5^(MHz), tan(77^) 




0,0 

H 

Cm 

Cl 

C 2 

min(5) 

(MHz) 

5 qubit 

95.8°, 86.3° 

-56.36, 

3.59, 

0.048 

209.48, 

22.51, 

0.0039 

-43.50, 

0.70, 

0.0029 

-37.56, 

0.70, 

0.004 

6 

3 qubit 

90°, 90° 

-56.0, 0, 

0 

211.83, 

8.98, 

0.0015 



56 


Table 4 Optimal orientations for the ENDOR control scheme with fully-labeled (5 qubit) and ^^Cm- 
labeled (3 qubit) molecules, corresponding values of hyperfine coupling constants and B^, and the 
quantity characterizing forbidden transition rates tan( 77 ”) for nuclear spin n. In these orientations, the 
NMR transitions are well separated for high fidelity control; the minimum distances between two NMR 
transition frequencies are 1.8 MHz and 7.6 MHz for the 5 qubit and 3 qubit systems, respectively. Here, 
min(5) represents the minimum separation between allowed ESR transitions. 


5 Simulation of Heat-Bath Algorithmic Cooling 

The Partner Pairing Algorithm (PPA) ca is the optimal method for implementing 
HBAC when there is one reset qubit that thermalizes to the bath polarization £ 7 ^ much 
faster than the relaxation rate of the system qubits. In this scenario, polarizing one 
system qubit beyond can be achieved by repeatedly applying the entropy compres¬ 
sion and refresh steps. The compression is a permutation that rearranges the diagonal 
elements of the density matrix in non-increasing order so that the polarization of the 
(target) qubit increases while the reset qubit polarization decreases. During the re¬ 
fresh step, the reset qubit thermalizes to the bath temperature and the overall entropy 
of the system is reduced. As this procedure is repeated, the polarization of the target 
qubit, asymptotically approaches a threshold value: £th = if ^ 2 “^^, and 

polarization arbitrarily close to unity when £ij ^ where n is the number of sys¬ 
tem qubits used in the compression step imiia. 
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We perform simulations of 3-qubit and 5-qubit HBAC for cooling one target qubit 
below the electronic spin bath temperature using -labeled and 1 ^ 2 -labeled 
molecules. The HBAC protocol consists of two steps: (i) apply a swap gate between 
the electron and each nucleus, with electron refresh steps in between, cooling all nu¬ 
clei to the bath temperature; (ii) run multi-qubit compression to boost the polarization 
of the target spin. For the 3 qubit case, the protocol we simulate is equivalent to the 
PPA. For the 5 qubit simulation, designing the stepwise operations corresponding to 
the PPA is challenging since each compression gate depends on the input state of 
the gate and is generally a different gate for each compression step. Moreover, the 
PPA compression gates may involve complicated pulse sequences. Thus, the 5-qubit 
protocol simulated here is designed to be analogous to the 3-qubit sequence, i.e. a po¬ 
larization transfer from the electron to all nuclei followed by a 5-qubit compression 
step. This same sequence can be repeated until asymptotic polarization is reached, 
although here we only simulated one round to demonstrate the feasibility of cooling 
beyond the Shannon bound in this 5-qubit system. 

Two control methods are tested for the HBAC simulations: AHC and ENDOR. 
The electronic T\ and T 2 processes are modelled as a Markovian dynamical map, and 
simulated by solving a master equation of the Lindblad form ||43l |44]|. Inhomoge¬ 
neous line broadening of the electron spin resonances, characterized by is taken 
into account by averaging the simulation over a set of spin Hamiltonians in which the 
magnitude of electron Zeeman energy is a Lorentzian-distributed random variable. 
We use experimentally measured electron 7i and T 2 values to determine the Lindblad 
operators, and the measured ESR line width to determine the T 2 Hamiltonian distri¬ 
bution. The electron relaxation parameters are 7i = 27 /is, 72 = 5 /is, and = 28 ns 
at room temperature. In the following sections, we describe the simulations in detail 
and report the results. 


5.1 Simulation of 3-qubit HBAC 

The quantum circuit for 3-qubit HBAC is shown in the top part of Fig. Here, the 
electron is chosen as the bath qubit whose polarization is 8 x 10“"^ at room tem¬ 
perature and at = 27r • 10 GHz/yg « 3568 G. Under these conditions, ^H and ^^C 
equilibrium polarizations are about 660 and 2620 times smaller than the electronic 
polarization. For the 3-qubit PPA with completely mixed initial states, the first two 
entropy compression steps are polarization transfers (swaps) from the reset qubit to 
the (target) and the 2^^ qubits. After the 3^^ gate operation (3-qubit compression), 
the target qubit polarization is higher than the bath polarization. These steps complete 
one round of the 3-qubit PPA. In theory, the first round of the 3-qubit PPA boosts the 
polarization of one nuclear spin to 1.5e^ — 0.5e^ 1.5e^ for small By repeated 

application of these steps, the nuclear spin is asymptotically polarized to 2e^. Note 
that starting from the 2^^ round, only one polarization swap gate is used since the 
target nucleus is already polarized above the bath polarization. 
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5.1.1 Anisotropic hyperfine control 

The crystal orientation used for the 3-qubit experiment using the AHC method is 
shown in Table In this orientation, the forbidden transition rate of in the pres¬ 
ence of a resonant microwave field is weaker than that of the a-proton. Consequently, 
the Cm polarization decays slower than that of proton during electron reset. Therefore 
Cm is chosen as the target qubit for cooling. The electron reset is done by waiting for 
4.2 X in order to bring the electron polarization to about 98.5% its thermal polar¬ 
ization. A longer waiting time, e.g. 5 brings the electron polarization to above 
99%, however, we find that the wait time of 4.2 x is optimal since nuclear spin po¬ 
larizations decay during electron reset (due to the anisotropic hyperfine interaction). 
The microwave field swap and compression gates are designed using the gradient as¬ 
cent pulse engineering (GRAPE) algorithm 1271 . The GRAPE pulse lengths are 840 
ns for the swap gates and 900 ns for the 3-qubit compression gate. The pulses are 
optimized over the electron Zeeman Hamiltonian distribution in order to be robust to 
= 28 ns. The design fidelities averaged over this distribution are 99%, 98% and 
96% for the e-H swap, e-C^ swap, and 3-qubit compression, respectively. The bottom 
part of Eig. [^illustrates the implementation of controls necessary for the first round 
of 3-qubit HBAC using AHC. 
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Fig. 6 Quantum circuit of 3-qubit HBAC using the AHC method. The electron spin is the reset qubit, and 
it is refreshed by waiting for 4.2 x Tf. In the optimal crystal orientation for AHC, the forbidden transition 
strength of Cm is weaker than that of the proton. Consequently, the Cm polarization decays more slowly 
during electron reset and Cm is chosen as the target qubit. Two rounds of the algorithm are shown in the 
top panel. The 3-qubit PPA iteratively applies polarization transfer (shaded in blue) and the compression 
(shaded in yellow). All gates are designed in the microwave domain using the GRAPE algorithm. The 
GRAPE pulse lengths are 840 ns for the swap gates and 900 ns for the compression gate. A schematic of 
the AHC sequence for this circuit is shown in the lower part. 


Eig.[7]shows the polarization of Cm at the end of each round of 3-qubit HBAC. 
The red curve is obtained when the simulation incorporates the experimentally deter¬ 
mined, room temperature electron relaxation effects. The Cm polarization exceeds the 
bath polarization after the first round of HBAC, and increases further asymptotically 
as the HBAC is repeated. Nevertheless, the polarization enhancement is below the 
theoretically calculated value; after 9 HBAC rounds, the polarization of Cm is about 
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Fig. 7 Simulation results for 3-qubit HBAC using AHC. The plot shows the ratio of the Cm polarization 
(Ec) and the electron bath polarization (e^) at the end of each HBAC round, up to 9 rounds. The black 
dashed curve is the theoretical (ideal) value. The red curve is obtained by incorporating all experimentally 
determined (room temperature) relaxation parameters. The blue curve is obtained by allowing the of 
the electron to equal T^. The green curve is obtained by allowing to be 20 /is, 4 times longer than the 
measured value, in order to test the consequence of a longer T 2 . 


76% of the theoretical value. The largest contribution to the error is the loss of nuclear 
polarization due to the electron T\ process during the application of GRAPE pulses 
and the electron refresh steps. In the absence of decoherence, the GRAPE pulses 
transfer 99.6% and 99.8% of the electron polarization to Cm and H, respectively. The 
compression gate polarizes the carbon to 1 A9£ij at the end of the first round. How¬ 
ever, when the T\ of the electron is introduced while T 2 and T 2 are still assumed to 
be infinite, the carbon and hydrogen polarizations prior to the compression step are 
reduced to 94.4% and 94.3% of the electron polarization, and the compression yields 
the carbon polarization of 1.39e/,. Another source of the error is the finite ratio of 
the electron T 2 to the pulse duration; T 2 is only about 5 to 6 times longer than the 
GRAPE pulses. One can imagine another type of 3-qubit molecule whose electron T 2 
is longer, for instance, 20 /is. Then the polarization of the carbon after 7 rounds of 
HBAC is 1.66£h, which is 83% of the theoretical value (i.e. the green curve in Eig.[^. 
Einally, we consider a scenario in which the ESR linewidth is much narrower and T 2 
limited, i.e. T 2 =T 2 = 5 /is. This result is indicated in blue in Eig.|^ Although there 
is a slight improvement by having a longer the experimental value of T 2 does 
not pose a significant problem since the GRAPE pulses are designed to be robust to 
inhomogeneous line broadening. 

5.7.2 ENDOR control 

In a pulsed ENDOR control scheme, the crystal orientation is chosen as shown in 
Table In this orientation, the forbidden transition rates are both minimized, and 
it turns out the rate for H is weaker than that of C^. Thus, H is chosen as the tar¬ 
get qubit for cooling. The electron reset is done by waiting for 5 x 7i in order to 
bring the electron polarization to 99.3% of the thermal polarization. The loss of 
nuclear polarization during reset and control operations is negligible since the for¬ 
bidden transition rates are both very weak. As shown in the bottom part of Eig. [8] 
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the swap and compression gates can be decomposed into controlled-not (CNOT) 
gates and a toffoli gate. These operations can be realized by selective microwave 
and RF ;r-pulses. For example, the CNOT gate that flips the H spin if the elec¬ 
tron is ‘spin down’ can be realized by RF pulses at frequencies that correspond to 
It//tele) ^ li/ztcie) and It/zicie) ^ \lHlcle) transitions. In the compression step, 
a Toffoli gate that flips H spin if both and electron are ‘spin down’ is required. 
However, the Toffoli gate on the proton cannot be realized by this method because the 
WH^cie) ^ transition frequency is identical to the |t//ic4'e) ^ li/rlcle) 

frequency. On the other hand, a Toffoli gate that flips the electron spin when both 
H and Cm are ‘spin down’ can be realized by applying a microwave pulse at the 
frequency of the \iHicte) ^ \iHicie) transition, which is distinct from all other 
allowed ESR transition frequencies. Therefore, we modify the quantum circuit to ex¬ 
tract the entropy from the electron during the compression, and use an additional 
swap gate to transfer the final polarization to the proton (see Fig.[^. 

This quantum circuit enables one to repeatedly apply HBAC for pumping the H 
polarization. The compression step consists of four CNOT gates and a Toffoli gate, 
as shown in Fig.[^ Since the goal of the compression is only to extract entropy from 
the electron, the last two CNOT gates targetting the nuclear spins (shown in the red 
dashed box in Fig. are not necessary. Hence, the compression is reduced to two 
CNOT gates and a Toffoli gate as shown in the middle part of Fig. ^ The rectangular¬ 
shaped microwave and RF pulses shown at the bottom of Fig. [8 represent selective 
;r-pulses, and the different shading of these pulses illustrate different frequencies. 


□= 


I = Polarization transfer = Compression 
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RF 
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Fig. 8 Quantum circuit of the 3-qubit HBAC using ENDOR control. The electron is used as the reset qubit, 
and it is refreshed by waiting for 5 x Ti. In the crystal orientation optimized for ENDOR, the forbidden 
transition rates of H are weaker than those of Cm, so H is chosen as the target qubit. The 3-qubit PPA 
iteratively applies polarization transfer (shaded in blue) and compression (shaded in yellow). All gates 
are realized using selective microwave and RF TT-pulses. The compression step can only be implemented 
here to boost the electron polarization (see text for explanation), so a swap gate between the electron and 
proton is used to store the boosted polarization on the proton. The CNOT gates in the red dashed box are 
not necessary and are left out of the simulated implementation. The quantum circuit in the middle shows 
the gate decomposition of the swap and the compression into CNOT and Toffoli gates. The schematic of 
the microwave and RF pulse sequence at the bottom shows that the CNOT and Toffoli gates can be realized 
by transition-selective pulses. Shading illustrates different pulse frequencies. 
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The maximum power for microwave pulses is chosen such that the Rabi frequency 
is 25 MHz. In practice, the resonant microwave cavity has a finite bandwidth. This 
means that more input power is needed in order to drive transitions whose frequency 
is offset with respect to the cavity resonance frequency. Since the input power is lim¬ 
ited in any real experimental setup, the pulse length must increase in order to apply 
the ;r-pulse at the offset frequency. Fig. shows a simulated ESR spectrum centered 
at the resonator resonance frequency of 10 GHz (red), and the voltage transfer func¬ 
tion for a resonator with quality factor (Q) of 100 (blue). If the pulse is applied at 
a frequency at which the corresponding value of the transfer function is v, then the 
pulse length must increase by a factor of 1/v to compensate for the loss in transmitted 
power. This finite bandwidth effect must be taken into account in simulations because 
the algorithm cannot be successful if microwave pulse durations become comparable 
to the electron T 2 . In the present simulations, the resonator quality factor Q is set to 
100 and microwave pulse lengths are adjusted accordingly. The RF ;r-pulses are re¬ 
alized using 15 /is (H) and 60 /is (C^) square pulses, refiecting typical RF amplifier 
output power levels and assuming an untuned (broadband) RF circuit. 



9.8 9.9 10 10.1 10.2 

Frequency (GHz) 


Fig. 9 Simulated field-sweep spectra for the ^^Cm-labeled malonic acid in the orientation given in Table|^ 
(red), and voltage transfer function for a microwave cavity with g = 100 (blue). The y-axis represents the 
scaling factor of square root of the microwave power in the cavity a function of the microwave frequency. 
Given a fixed maximum available microwave power, pulse durations for offset pulses must be increased 
relative to cavity-resonant pulses. 


Using the room-temperature electron Ti, the HBAC algorithm will not be suc¬ 
cessful because RF pulse lengths are similar to Tj. This can be solved by exploiting 
longer Tf values at lower temperatures. In unlabelled, irradiated malonic acid, we 
have found experimentally a Ti that grows roughly exponentially with temperature, 
e.g. with values of 29 /is, 2.6 ms and 11 ms at room temperature, 43 K, and 22 K, 
respectively. For simulating 3-qubit HBAC with ENDOR, we choose Tj = 2.6 ms 
(ignoring the fact that a structural phase transition occurs at 47 K, probably compli¬ 
cating the ESR spectrum). The selective microwave pulses must be designed with 
a care. While each pulse selectively excites certain transitions, they must also be 
broad enough to cover the ESR linewidth of about 12 MHz. First, we simulate 50 
ns gaussian-shaped pulses with the full width at half maximum of 20 ns in order 
to excite the entire ESR linewidth and remain selective on the particular transition. 
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The result is shown in Fig. in red. The main source of the error is the inability of 
the Gaussian pulses to uniformly rotate all spins across the ESR linewidth; the pulse 
bandwidth must be close to 50 MHz in order to fully excite the entire ESR line width, 
but the minimum distance between two ESR allowed transition frequencies is about 
56 MHz. Thus, even in the absence of T\ and T 2 effects, each swap gate loses 8% 
polarization. For the improved microwave control, selective ;r-pulses are engineered 
using the GRAPE algorithm and can be made robust to T 2 =28 ns which corresponds 
to the ESR line width of 12 MHz. Using this method, swap gates can transfer 98% 
of the polarization from the electron to the target nuclear spin. The simulation results 
for HBAC using GRAPE microwave pulses and rectangular RF pulses are shown in 
blue in Fig.p^ 



Fig. 10 Simulation results of the 3-qubit HBAC using ENDOR control. The plot shows the ratio of H 
polarization to the electron thermal polarization after each round of the HBAC, up to 8 rounds. The black 
dashed curve is the theoretical value. The selective microwave TT-pulses are designed in two different ways: 
(i) 50 ns gaussian-shaped MW pulses with 20 ns full width at half maximum (red), (ii) and GRAPE pulses 
(blue). In both cases, RE transitions are applied using square transition-selective pulses with 15 ns and 60 
ns pulse lengths for H and Cm spins, respectively. 


5.2 Simulating 5-qubit HBAC 


As discussed in Section 4.1 due to weak forbidden transition rates of Ci 9 and fre 


quency overlap among ESR transitions, the AHC scheme cannot be implemented in 
the 5-qubit sample. Therefore, we focus on ENDOR control techniques in the simu¬ 
lation of 5-qubit HBAC. Here, we simulate one round of the 5-qubit HBAC as shown 
in Fig. m to demonstrate that cooling beyond the Shannon bound is experimentally 
feasible in the 5-qubit molecule. The quantum circuit we employ for 5-qubit HBAC 
is shown in Fig.[TT] The electron is the target qubit, but the circuit can be easily mod¬ 
ified to cool a nuclear spin by adding a swap gate after the compression step. The¬ 
oretically, by applying one round of the quantum circuit, the target qubit is cooled 
to £^(15 — 10e| -j- 3e^)/8 ^ 1.875e^ for small The orientation of the magnetic 
field used in the simulation is shown in Table and the electron polarization and 
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relaxation parameters are the same as used in the 3-qubit ENDOR HBAC simulation. 
The forbidden transition rates of the carbons are weaker than those of the proton, so 
the polarization transfer to H is done just before the compression gate, as shown in 
Fig. [TT] Similar to the 3-qubit HBAC ENDOR circuit, the last CNOT gates in the 
compression step (inside the red dashed box) are not necessary since the electron is 
the target spin. 


e 

Cm 

Cl 

C2 
H 

Fig. 11 Quantum circuit for 5-qubit HBAC with the electron as the target spin. The polarization can be 
transferred to a nuclear spin with an additional swap gate. The last CNOT gates inside the red dashed 
box of the compression step are unnecessary, and are left out of the simulation. The open circles in the 
controlled gates indicate that the target spin is flipped if the control qubit is in the ‘spin up’ state. One 
round of the algorithm shown. 
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The usual method for a swap gate between the electron and one of the carboxyl 
carbons requires selective microwave control of Ci and C 2 . However, the hyperfine 
tensors of Ci and C 2 are very similar and the two spins cannot be separately addressed 
by microwave pulses. However, they can be separately addressed with RE pulses 
(see Section [ 4 ^. Therefore, the polarization transfer step for Ci^2 is modified as 
shown in Fig.|12[ and theoretically polarizes both Ci and C 2 to the electron thermal 
polarization. 





Fig. 12 Circuit for polarization transfer between the electron and Ci^2 • The circuit brings Ci and C 2 to the 
bath temperature without selective microwave control, but requires selective RF pulses. The open circles 
in the controlled gates indicate that the target spin is flipped if the control qubit is in the ‘spin up’ state. 


Note that the circuit shown in Fig. [TT]is not the optimal cooling algorithm, i.e. 
the PPA. Fig. [^compares the PPA and our cooling algorithm by showing the target 
qubit polarization £t compared to the bath polarization as a function of HBAC steps. 
Here, each step consists of one refresh operation and one gate on the system qubits (in 
contrast to the ‘round’ used previously). In the low bath polarization regime, the PPA 
asymptotically increases the target qubit polarization to 8 times the bath polarization, 
while repeated application of the algorithm shown in Fig. [TT] yields an asymptotic 
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enhancement of 4. While it is possible to find sequences corresponding to the PPA, it 
is an open question whether such sequences could be practically implemented in this 
system. For the PPA, the gate decomposition of each system qubit operation depends 
on the input state, and can result in a complicated pulse sequence. An advantage of 
the algorithm implemented in our simulation is that the gate decomposition of each 
step is relatively simple and the asymptotic limit can be reached by simply repeating 
the quantum circuit. 



Fig. 13 Theoretical target qubit polarization normalized by the bath polarization, as a function of the 
number of cooling steps, to show the difference between the 5-qubit PPA and the algorithm in Fig. m 
Each step consists of one refresh operation and one gate on the system qubits. The red curve is obtained by 
the PPA while the blue curve is obtained by repeated application of the quantum circuit shown in Fig. m 
In the low polarization regime, the PPA allows the target qubit polarization to asymptotically approach 8 
times the bath polarization, and 4 times the bath polarization with our algorithm. 


As seen already in the 3-qubit HBAC simulation, designing microwave pulses 
that uniformly rotate all spin across the ESR line width while remaining transition- 
selective is a challenge. For the 5-qubit system, this problem is exacerbated since the 
ESR transitions are more closely spaced in frequency. In order to obtain the maximum 
polarization enhancement, GRAPE optimal control must again be used for designing 
the selective microwave pulses. GRAPE pulses that are robust to the experimentally 
determined value of have been designed with at least 0.95 state fidelity and with 
200 ns pulse lengths for the mw pulses applied to Cm and C 2 , 900 ns for Ci and ^H, 
and 100 ns for the 5-qubit compression. RE pulse lengths are chosen as 5 /is and 20 
/is for and respectively. Including the effects of the electron 7i, T 2 , and 
the polarization enhancement after one round of 5-qubit HBAC is 1.67e^. A major 
contribution to the error is the broad ESR line width compared to the spacing between 
ESR resonances, which makes realizing high fidelity transition-selective rotations on 
a time scale short compared to T 2 difficult. Despite this, the simulation shows that a 
single round of 5-qubit HBAC with a realistic control sequence and relaxation yields a 
polarization enhancement beyond the Shannon bound, and also beyond the capability 
of a single round of 3-qubit HBAC. Another simulation is performed assuming a 
much narrower ESR line width, i.e. T 2 = T 2 . In this case, all microwave pulses can 
be designed with 99% unitary fidelity at 100 ns pulse length, and all nuclei reach 
at least 97% of the electron polarization during the polarization transfer steps. The 
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polarization improvement after one round of 5-qubit HBAC is 1.79e^. This is 95% of 
the polarization improvement predicted by the theory, similar to the result of the 3- 
qubit HBAC simulation using ENDOR control with GRAPE microwave pulses. We 
expect the single qubit polarization to reach above 3.6e/,, 90% of the theoretical value, 
as the algorithm is repeated, similar to the 3-qubit case. The result shows that for 
HBAC with many qubits it is crucial to have a sharp ESR linewidth, which motivates a 
search for electron-nuclear spin systems having narrower ESR lines than the malonyl 
radical. 

The 5-qubit HBAC results are summarized in Table 




Simulation 


Theory 

II 

Ty = 28 ns 

etl£b 

1.875 

1.79 

1.67 


Table 5 Polarization of the target qubit compared to the bath qubit after one round of 5-qubit HBAC. The 
quantum circuit for one round of 5-qubit HBAC is shown in Fig. Results are shown for two different 
values of T^. The electron spin lattice relaxation time is Ti =2.6 ms in this simulation. 


Einally, we note two experimental caveats that are not taken into account in our 
simulations, but do not invalidate the main results. In section we accounted for 
a type of cross-relaxation in which the nuclear polarization decays due to a com¬ 
bination of the electronic spin-lattice relaxation process and the anisotropic term in 
the spin Hamiltonian of the form SJx- There is an additional cross-relaxation mech¬ 
anism involving noise acting directly on the nuclear operators, which we have not 
included in the Lindblad master equation. In future work, we plan to experimentally 
measure the cross-relaxation rates as a function of orientation and temperature to 
determine these additional contributions. If the additional mechanism is dominant, 
going to lower temperatures may be required in order to achieve sufficiently long 
nuclear Ti timescales. Secondly, a well known structural phase transition takes place 
in malonic acid at 46 K |[34l - [36ll . which we have not considered here. Below 46 K, 
the PI crystal symmetry is broken and there are two magnetically distinguishable 
molecules per unit cell. We have observed certain ENDOR transitions to split into 
two below 46 K, consistent with this phase change. However, experiments may be 
designed using phase-cycling techniques in order to cancel signal contributions from 
one of the molecules in the unit cell, so this does not in principle prevent the proposed 
AC experiments from being carried out. 


6 Conclusion 

HBAC is an open-system cooling method that allows at least one qubit to be polarized 
beyond the Shannon bound, starting from a set of initially mixed qubits and exploit¬ 
ing a controlled interaction with a heat-bath. This is a promising tool for dynamically 
preparing ancilla qubits to a sufficient level of purity for quantum error correction. In 
this work, we accurately determined the full spin Hamiltonian of a five qubit electron- 
nuclear hyperfine coupled system in single crystal, irradiated malonic acid. Using the 













Hyperfine spin qubits in irradiated malonic acid: heat-bath algorithmic cooling 


23 


hyperfine tensors, we determined the optimal magnetic field orientation for achieving 
high fidelity control using two control methods: anisotropic hyperfine control (AHC) 
and pulsed ENDOR techniques. Computer simulations were carried out using real¬ 
istic experimental conditions and including the relevant electron spin relaxation pa¬ 
rameters, demonstrating that the realization of 3-qubit and 5-qubit HBAC is feasible 
in this system. Using the 3-qubit molecule with AHC, the polarization of a nuclear 
spin is predicted to increase above 1.5e/, after 9 rounds of the cooling algorithm. The 
ENDOR simulation assumed a lower temperature in which the electron T\ =2.6 ms, 
in order to have the Ti sufficiently long compared to RE pulse lengths. Using GRAPE 
for selective microwave pulse design, the polarization of a nuclear spin is predicted 
to increase above 1.8e^ after 9 rounds of the cooling algorithm. There is a tradeoff 
here between experimental simplicity and the achievable fidelity; the AHC experi¬ 
mental setup is simpler and can be performed at room temperature, but in general 
yields lower fidelities than the ENDOR approach at low temperature. 

The carboxyl carbons in the 5-qubit sample have very weak forbidden transition 
rates under an applied microwave field, which prevents them from being practically 
controllable using AHC. Hence, for the 5-qubit HBAC simulation, we focused on the 
ENDOR control scheme. The 5-qubit simulation was again carried out at low tem¬ 
perature for a longer T \, and the selective microwave pulses were designed using the 
GRAPE method. After one round of 5-qubit HBAC, the simulation predicts that the 
target qubit (the electron) would reach a polarization of 1.67The major obstacle to 
reaching the ideal value of 1.875e^ is the small spacing between certain allowed ESR 
transitions, similar in order to the ESR line width. Another simulation was performed 
assuming = T 2 to predict the outcome of the algorithm given a molecule that is 
similar to the malonic radical but with much sharper ESR transitions. Here, a target 
qubit polarization of 1.79e^ is obtained after one round of the algorithm. 

We conclude that the experimental demonstration of 3- and 5-qubit algorithmic 
cooling beyond the Shannon bound is feasible in the isotopically labelled malonyl 
radical. The 5-qubit system can yield a larger polarization enhancement compared 
to the 3-qubit system, as expected. The experimental value of T 2 is found to be a 
critical factor that limits the fidelity of gate operations, and therefore the achievable 
polarization under HBAC. Nearly ideal results are obtainable when T 2 = T 2 . 
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